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ABSTRACT 

We present a systematic analysis of archival HST WFPC2 "Association" data sets that correlate 
with the Chandra positions of a set of 44 ultra-luminous X-ray sources (ULXs) of nearby galaxies. The 
main motivation is to address the nature of ULXs by searching for optical counterparts. Sixteen of the 
ULXs are found in early- type galaxies (RC3 Hubble type < 3). We have improved the Chandra-HST 
relative astrometry whenever possible, resulting in errors circles of 0.3-1.7" in size. Disparate numbers 
of potential ULX counterparts are found, and in some cases none are found. The lack of or low number 
of counterparts in some cases may be due to insufficient depth in the WFPC2 images. Particularly in 
late-type galaxies, the HST image in the ULX region was often complex or crowded, requiring source 
detection to be performed manually. We therefore address various scenarios for the nature of the ULX 
since it is not known which, if any, of the sources found are true counterparts. The optical luminosities 
of the sources are typically in the range 10'^~^Lq, with (effective) V magnitudes typically in the range 
22-24. In several cases color information is available, with the colors roughly tending to be more red 
in early-type galaxies. This suggests that, in general, the (potential) counterparts found in early-type 
galaxies are likely to be older stellar populations, and are probably globular clusters. Several early- 
type galaxy counterparts have blue colors, which may be due to younger stellar populations in the 
host galaxies, however these could also be background sources. In spiral galaxies the sources may also 
be due to localized structure in the disks rather than bound stellar systems. Alternatively some of the 
counterparts in late-type galaxies may be isolated supergiant stars. The observed X-ray/optical flux 
ratio is diluted by the optical emission of the cluster in cases where the system is an X-ray binary in a 
cluster, particularly in the case of a low-mass X-ray binaries in old cluster. If any of the counterparts 
are bound systems with ~ 10*~^ stars and are the true counterparts to the ULX sources, then the 
X-ray luminosities of the ULX are generally well below the Eddington limit for a black hole with mass 
^ 0.1% of the cluster mass. Finally, we find that the optical flux of the counterparts is consistent with 
being dominated by emission from an accretion disk around an intermediate-mass black hole if the 
black hole happens to have a mass > 1O^M0 and is accreting at close to the Eddington rate, unless 
the accretion disk is irradiated (which would result in high optical disk luminosities at lower black 
hole masses). 

Subject headings: catalogs — X-rays: binaries — X-rays: galaxies 



1. INTRODUCTION 

In the early 1980s, surveys of normal galaxies with the 
Einstein X-ray satellite revealed intermediate-luminosity 
(Lx ~ lO'^^ — 10^" erg s~^) X-ray sources which were 
seemingly located in the centers of spiral galaxies (Fab- 
biano 1989). This was very interesting, since Seyfert nu- 
clei (active galactic nuclei [AGNs] in nearby spirals) are 
typically much more luminous (Lx ^ 10^^ — lO*** erg s""'^) 
and black hole X-ray binaries (BH XRBs) are much less 

^ Based on observations made with the NASA/ESA Hubble 
Space Telescope, obtained from the Data Archive at the Space 
Telescope Science Institute, which is operated by the Association 
of Universities for Research in Astronomy, Inc., under NASA con- 
tract NAS 5-26555. This project is associated with Archival pro- 
posal #9545. 



luminous (typically Lx ^ 10"^^ erg s^^). It was not 
clear whether these intriguing sources were underlumi- 
nous accreting supermassive BHs, overluminous XRBs 
located near the galactic nucleus, or an entirely new 
type of astrophysical object altogether (e.g., Colbert et 
al. 1995). In the 1990s, ROSAT High Resolution Im- 
ager (HRI) observations showed that these Ultralumi- 
nous X-ray objects (ULXs; also known as Intermediate- 
luminosity X-ray Objects, or IXOs) are compact X-ray 
sources, and are quite common in the local Universe 
(Colbert & Mushotzky 1999; Roberts & Warwick 2000). 
Many ULXs were found to be offset from the optical nu- 
cleus (e.g., Colbert & Mushotzky 1999). More complete 
surveys from ROSAT (e.g., Roberts & Warwick 2000; 
Colbert & Ptak 2002) and Chandra imaging data place 
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their frequency of occurrence at about one in every five 
galaxies (Ptak & Colbert 2004). Most of the ULXs found 
in these surveys are clearly significantly displaced from 
the galactic nuclei, implying that their masses must be 
<10^ Mq since dynamical friction would otherwise cause 
them to sink to the center of the galaxy in a Hubble time 
(Tremaine, Ostriker & Spitzer 1975). 

ULXs are generally unresolved at the high spatial res- 
olution (~ 0.5") of Chandra, and many show variability 
(Colbert & Ptak 2002; Fabbiano et al. 2003), which rules 
out the hypothesis that ULXs are closely-spaced aggre- 
gates of lower-luminosity X-ray sources. Bondi accretion 
from a dense ISM may not be sufficient to explain the 
observed luminosities (King et al. 2001), although an in- 
termediate mass black hole (M > 300Mq) in a dense 
molecular cloud may achieve luminosities in excess of 
]^q39 ^,j,gy y-1 (Ki-Qii]j 2004) . They arc therefore most 

likely powered by accretion in binary systems. This is 
supported by evidence for periodicity in a few individ- 
ual objects (Sugiho et al. 2001; Liu et al. 2002b). If one 
assumes that ULXs have Lx below the Eddington lumi- 
nosity and are not beamed, then the mass of the central 
object is required to be in the range ^ 15 — 500Mq. 
This is too massive to be a neiitron star or a black 
hole formed through normal stellar evolution, thus we 
would have to be dealing with intermediate-mass black 
holes (IMBHs). This interpretation is quite fascinating, 
as reviewed by van der Marel (2004) and Miller & Col- 
bert (2004). For some individual ULXs it has been ar- 
gued that the presence of an IMBH is plausible (Porte- 
gies Zwart et al. 2004), but in general it poses many 
challenges from a theoretical point of view (see King 
et al. 2001; King 2002). Alternatively, ULXs may in- 
stead be stellar-mass black hole systems that arc emit- 
ting anisotropically with either mild (King et al. 2001) or 
relativistic beaming (Kording, Falcke & Markoff 2002), 
or are actually emitting at super-Eddington rates (Begel- 
man 2002; Grimm, Gilfanov & Sunyaev 2002). For ex- 
ample, the galactic microquasar GRS 1915 has a known 
black hole mass of 14M0 and is often observed with X- 
ray luminosities exceeding 10^^ ergs s~^, qualifying it as 
an ULX. In this case the high luminosity is due to both 
relativistic beaming and super-Eddington emission (King 
2002; Mirabel & Rodriguez 1999). 

Phenomenologically there appear to be several breeds 
of ULXs. They may be preferentially found in star- 
forming galaxies (Kilgard et al. 2002; Grimm et al. 2003), 
with starburst galaxies such as the Antennae (Zezas et 
al. 2002; Zezas & Fabbiano 2002) and the Cartwheel 
(Gao et al. 2003) having large numbers of ULXs. This 
suggests an association between the incidence of ULXs 
and star formation (although Ptak & Colbert [2004] have 
found that this is not necessarily the case in general) . In- 
deed, optical counterparts reported for ULXs suggest a 
young star cluster in one case (Goad et al. 2002) and a 
single 0-star in another case (Liu, Bregman & Seitzer 
2002a). Terashima, Inoue, & Wilson (2006) discuss a 
study of optical counterparts to ULXs in M51 and vary- 
ingly find no, one or several candidate counterparts. The 
colors and absolute magnitudes of counterparts in several 
cases (where there were only one or two counterparts) 
were consistent with stars of mass 7-9 Mq and 10-15 
Mq, implying that these ULXs may be high-mass X-ray 
binaries. A radio counterpart was found for a ULX in 



the dwarf galaxy NGC 5408 (Kaarct et al. 2003), with 
two possible HST counterparts in the X-ray error cir- 
cle. The brighter of these counterparts has My ~ —6.2, 
consistent with a supergiant O star. 

ULXs have also been found in elliptical galaxies and 
galactic bulges (Colbert & Ptak 2002; Angelini et 
al. 2001; Wu et al. 2002; Jeltema et al. 2003). This 
suggests that some ULXs are associated with old sys- 
tems such as globular clusters. However, ULXs with 
Lx > 2 X 10^^ ergs s~^ found in ellipticals are often 
not associated with globular clusters and are observed in 
numbers consistent with the expected number of back- 
ground sources (Irwin et al. 2004). A study of optical 
counterparts to four ULXs by Gutierrez (2006) found 
that three are background AGN and one is likely to be 
foreground star. Three of these ULXs were found in 
early-type galaxies, again showing a tendency for ULXs 
associated with early-type galaxies to be background 
AGN. 

X-ray spectroscopy can give clues about accretion disks 
and coronae that may be present in ULXs (Colbert & 
Mushotzky 1999; Makishima et al. 2000; Roberts et 
al. 2001; Strickland ct al. 2001; Zezas et al. 2002; Fos- 
chini ct al. 2002; Miller ct al. 2003). However, there are 
important questions about the nature of ULXs that can 
only be answered using the multi-wavelength properties 
of as large a sample of ULXs as possible. In the Anten- 
nae, ULXs are often observed close to star clusters, but 
not coincident with them. This may suggest a scenario 
in which ULXs are X-ray binaries that can be ejected out 
of clusters through recoil (Portegies Zwart & McMillan 
2000). If this is generally true then it would preclude 
an IMBH because the mass would be too large for the 
binary to be ejected (Miller & Hamilton 2002). In cases 
where individual counterparts are not detected learning 
about the environment of the ULXs can also be instruc- 
tive. For example, it has been suggested that ULXs in 
late-type galaxies arc preferentially found in or near HII 
regions (PakuU & Mirioni 2002). 

We have established a pipeline to continuously analyze 
Chandra data from galaxies as it becomes publicly avail- 
able, with emphasis on producing a catalog of ULXs. The 
present paper presents a study of archival Hubble Space 
Telescope (HST) images obtained with the Second Wide 
Field and Planetary Camera (WFPC2) of fields that eon- 
tain ULXs in our Chandra catalog. The goal is to search 
for optical counterparts to the ULXs, and to character- 
ize their environment. While a few individual optical 
ULX coimtcrparts have been reported by other groups, 
the present work represents the first large-scale Chan- 
dra/HST ULX cross-comparison and cataloging effort. 
The use of HST data is critical, given that ground-based 
observations do not have the spatial resolution to resolve 
individual stars (unless the star is unusually isolated) 
and often cannot reach the apparent magnitudes of even 
supergiant stars except for the nearest galaxies. The re- 
sults provide much improved statistics with fewer biases 
than previous work, as well as a catalog of sources for 
use in follow-up studies. 

The paper is structured as follows. Section 2 discusses 
the sample selection and the Chandra and HST data that 
we have used for the study. Section 3 discusses the astro- 
metric registration of the Chandra and HST data. Sec- 
tion 5 discusses the results from photometric analysis 
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on the potential counterparts that were located. Sec- 
tion 6 discusses the main results of the study, and the 
implications for our understanding of ULXs. Section 7 
summarizes the main conclusions. 

2. SAMPLE SELECTION 
2.1. Catalog of Chandra ULXs 

We searched for Chandra X-ray sources within the R25 
radius of all galaxies in the Third Reference Catalog of 
Bright Galaxies (RC3; de Vaucouleurs et al. 1991) with 
recessional velocities cz < 5000 km s~^. Here R25 is the 
25 mag arcsec"^ isophotal diameter, as listed in RC3. 
An X-ray point source database was created with the 
XAssist (Ptak & Griffith 2003) Chandra ACIS pipeline, 
using all public ACIS data as of 17 July 2002. ULXs 
were selected by requiring a 2—10 keV X-ray luminos- 
ity Lx > 10^^ ° erg s~^, assuming an absorbed power- 
law model with Galactic absorption (Dickey & Lockman 
1990) and photon index T — 1.7. For galaxies with 
RC3 recessional velocities cz < 1000 km s~^, distances 
were taken from TuUy (1988), if available, otherwise from 
Kraan-Kortewag & Tamman (1979). Distances to galax- 
ies with cz > 1000 km were computed using Ho = 75 
km Mpc^^ (for consistency with TuUy 1988). H I line 
velocities from RC3 were used when available, otherwise 
"optical" RC3 velocities were used. 

The preliminary source list thus obtained was tapered 
by merging duplicate sources and examining the indi- 
vidual X-ray images to verify that the source was not 
an X-ray instrumental artifact and that it was indeed 
point-like. We next cross-correlated the X-ray positions 
with known sources listed in the NASA Extragalactic 
Database (NED) to identify known AGNs, galaxy nuclei, 
and X-ray supernovae. Any objects thus identified were 
omitted, and questionable objects were flagged as such. 
The remaining list contained 125 ULXs, which were used 
as a starting point for searching the HST WFPC2 archive 
using the selection criteria discussed in Section 2.2. For 
those ULXs with WFPC2 data that passed our search 
criteria, we performed an astrometric registration as de- 
scribed in Section 3. With more accurate absolute po- 
sitions available, we reexamined all of the questionable 
objects by checking their X-ray position with respect to 
the galaxy nucleus in the WFPC2 images, and by reject- 
ing sources when spatially coincident with known AGNs 
or SNe. This ultimately yielded 44 ULXs in 25 different 
galaxies, that make up our sample for this paper. Table 1 
lists the basic properties of the sample. 

2.2. Selection of HST/WFPC2 Data 

For the HST parts of our project we chose to work with 
WFPC2 associations, instead of individual WFPC2 im- 
ages. Associations are logically related sets of individual 
WFPC2 images (e.g. the same field, filter, and orienta- 
tion) with accompanying information on positional off- 
sets between the images. The vast majority of WFPC2 
science data (93%) has been made into associations (Mi- 
col & Durand 2002). The WFPC2 Associations Science 
Products Pipeline^ (WASPP) produces data products 

^ WASPP was developed by the Canadian Astronomy Data Cen- 
tre (CADC) in eonjunction with the Space Telescope European 
Coordinating Facility (ST-ECF) and the Multimission Archive at 
STScI (MAST). 



from the images in an association. To this end it starts by 
applying the regular WFPC2 pipeline calibration steps 
(bias subtraction, dark subtraction, flat-fielding, etc.) to 
the individual images of the association. The reduced 
images are then registered and co-added with cosmic-ray 
rejection. If reference stars from the USNO-A2.0 Catalog 
(Monet et al. 1998) are present in the field, then these 
are identified by an automated algorithm and the final 
combined image is registered to their reference frame. 
For the present project we used the data made available 
in the first release (November 6, 2002) of the WASPP 
association data. This release contained data products 
for 66% of all WFPC2 associations. 

We searched the STScI HST Data Archive for WFPC2 
images with pointings that might include the ULXs 
that were identified with Chandra (Section 2.1). Each 
WFPC2 field was then visually inspected to ensure that 
the ULX position fell on one of the WFPC2 chips. Fields 
in which the ULX position was located off the edge of the 
WFPC2 detectors were eliminated. Of these images, 90% 
was part of an association. Each association typically 
contains 2 8 images. Of these associations, 74% had been 
processed by the WASPP at the time of the first data re- 
lease. This implies that two-thirds of all the images that 
we could have used for the present study were part of 
the first WASPP release. Since manual registration and 
co-addition of the remaining one-third would have been a 
substantial effort while providing only a modest increases 
in overall sample size, we decided to ignore these data for 
the present study. The final WFPC2 data sample thus 
obtained consists of 76 WFPC2 images that include the 
44 ULXs in our Chandra catalog. Each WFPC2 image 
is a fully calibrated association data product. For a de- 
tailed description of the WASPP data reduction steps we 
refer the reader to Schade et al. (2002). The WFPC2 as- 
sociation products used in our study are listed in Table 2, 
together with basic properties such as the filter used and 
the total exposure time. For some of our ULXs we have 
images in multiple filters, or even multiple images in the 
same filter (the WASPP pipeline does not combine im- 
ages taken at different epochs and telescope orientation, 
even when taken with the same filter). 

3. ASTROMETRIC REGISTRATION 

It is important for our study that the HST and Chan- 
dra data are astrometrically aligned as well as possible. 

The absolute astrometric accuracy of WFPC2 data is 
dominated by two components. The first is the ^ 0.5" 
RMS astrometric accuracy of the Guide Star Catalog, 
which is used for pointing and guiding. The second is the 
accuracy with which the relative positions of the WFPC2 
and the Fine Guidance Sensors in the focal plane are 
known. The latter effect dominates the overall error bud- 
get, although its actual size has varied over the WFPC2 
lifetime (Brammer et al. 2002). By comparing the posi- 
tions of USNO-A2.0 stars in large numbers of WFPC2 
images to their known astrometric coordinates, Schade et 
al. (2002) found an overall RMS astrometric accuracy for 
WFPC2 data of 1.6". By contrast, Chandra has nominal 
90% and 99% confidence astrometric accuracies of 0.6" 
and 0.8" (Aldcroft 2002). 

The WASPP association products have the advantage 
that they have significantly enhanced astrometric accu- 
racy as compared to the nominal WFPC2 accuracy. The 
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TABLE 1 
ULX Sample 



ULX 


RA 


DEC 


Galaxy 


T-type 


Dist. 


Chandra 


Lx 


ID 




(J2000) 




(RC3) 


(Mpc) 


ObsID 


(1039 erg s-i) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 



1 


0:52:02.72 


47:33:07.7 


NGC 278 


3.0 


11.8 


2056 


2.0 


2 


2:22:31.36 


42:20:24.4 


NGC 891 


3.0 


9.6 


794 


1.8 


3 


2:42:38.89 


-00:00:55.1 


NGC 1068 


3.0 


15.2 


344 


2.9 


4 


2:42:39.71 


-00:01:01.4 


NGC 1068 


3.0 


15.2 


344 


1.0 


5 


2:46:19.82 


-30:16:02.6 


NGC 1097 


3.0 


17.0 


1611 


2.9 














2339 


2.8 


6 


3:38:27.64 


-35:26:48.2 


NGC 1399 


-5.0 


19.3 


319 


1.0 


7 


3:38:31.81 


-35:26:03.8 


NGC 1399 


-5.0 


19.3 


319 


1.3 


8 


3:38:32.60 


-35:27:05.1 


NGC 1399 


-5.0 


19.3 


319 


1.6 


9 


4:56:52.21 


-4:52:17.9 


NGC 1700 


-5.0 


52.2 


2069 


2.3 


10 


4:56:56.03 


-4:51:59.5 


NGC 1700 


-5.0 


52.2 


2069 


3.0 


11 


4:57:01.82 


-4:51:15.9 


NGC 1700 


-5.0 


52.2 


2069 


6.4 


12 


9:55:46.44 


69:40:40.5 


NGC 3034 


90.0 


5.2 


379 


1.0 


13 


9:55:50.01 


69:40:46.0 


NGC 3034 


90.0 


5.2 


361 


2.1 














1302 


2.2 














378 


6.9 














379 


9.0 


14 


9:55:50.91 


69:40:46.6 


NGC 3034 


90.0 


5.2 


378 


2.8 


15 


10:27:52.54 


-43:53:50.2 


NGC 3256 


99.0 


37.1 


835 


1.3 


16 


10:27:55.15 


-43:54:47.2 


NGC 3256 


99.0 


37.1 


835 


2.3 


17 


10:36:42.66 


-27:31:40.8 


NGC 3311 


-4.0 


50.5 


2220 


4.5 


18 


10:36:42.75 


-27:31:43.5 


NGC 3311 


-4.0 


50.5 


2220 


9.2 


19 


11:20:15.76 


13:35:13.7 


NGC 3628 


3.0 


7.7 


2039 


1.2 














395 


2.3 


20 


11:20:20.90 


12:58:46.0 


NGC 3627 


3.0 


6.6 


394 


1.2 


21 


12:01:51.32 


-18:52:25.4 


NGC 4038/9 


9.0 


21.7 


315 


2.0 


22 


12:01:52.08 


-18:51:33.6 


NGC 4038 


9.0 


21.7 


315 


5.5 


23 


12:01:54.27 


-18:52:01.9 


NGC 4038/9 


9.0 


21.7 


315 


1.0 


24 


12:01:54.35 


-18:52:10.3 


NGC 4039 


9.0 


22.1 


315 


1.0 


25 


12:01:54.97 


-18:53:15.1 


NGC 4038/9 


9.0 


21.7 


315 


2.8 


26 


12:01:55.65 


-18:52:15.1 


NGC 4038/9 


9.0 


21.7 


315 


3.8 


27 


12:01:56.43 


-18:51:57.9 


NGC 4038/9 


9.0 


21.7 


315 


3.6 


28 


12:01:58.22 


-18:52:04.5 


NGC 4038/9 


9.0 


21.7 


315 


1.0 


29 


12:10:33.77 


30:23:57.9 


NGC 4150 


-2.0 


9.7 


1638 


1.1 


30 


12:30:43.26 


41:38:18.4 


NGC 4490 


7.0 


7.8 


1579 


1.1 


31 


12:35:58.56 


27:57:41.9 


NGC 4559 


6.0 


9.7 


2026 


2.2 














2027 


4.2 


32 


12:36:17.40 


25:58:55.5 


NGC 4565 


3.0 


16.4 


404 


9.8 


33 


12:40:00.35 


-11:37:24.0 


NGC 4594 


1.0 


14.5 


407 


1.4 


34 


12:50:25.68 


25:31:29.2 


NGC 4725 


2.0 


16.1 


409 


1.6 


35 


13:05:21.94 


-49:28:26.6 


NGC 4945 


6.0 


5.2 


864 


1.0 


36 


13:12:55.59 


-19:30:39.7 


NGC 5018 


-5.0 


37.3 


2070 


5.3 


37 


13:13:02.03 


-19:31:05.5 


NGC 5018 


-5.0 


37.3 


2070 


2.5 


38 


13:13:29.66 


36:35:23.1 


NGC 5033 


5.0 


18.7 


412 


2.2 


39 


13:13:29.46 


36:35:17.4 


NGC 5033 


5.0 


18.7 


412 


1.5 


40 


13:25:19.84 


-43:03:17.2 


NGC 5128 


-2.0 


4.9 


316 


1.0 


41 


14:13:10.03 


-65:20:45.1 


Circinus 


3.0 


3.7 


355 


1.6 


42 


14:13:12.24 


-65:20:14.0 


Circinus 


3.0 


3.7 


356 


1.8 














365 


4.9 


43 


17:56:01.59 


18:20:22.6 


NGC 6500 


1.7 


40.0 


416 


2.5 


44 


22:02:07.92 


-31:59:19.7 


NGC 7174 


1.8 


37.0 


905 


1.5 



Note. — The table lists basic properties for the ULXs in our sample, selected as described in Section 2.1. Column (1) lists the running ID number 
of the ULX used throughout this paper. Columns (2) and (3) list the X-ray source position from the Chandra ACIS data. Column (4) lists the 
host galaxy name, and column (5) lists the galaxy morphology T-type from the Third Reference Catalog of Bright Galaxies (RC3, de Vaucouleurs 
et al. 1991). Column (6) lists the galaxy distance, obtained as described in Section 2.1. Column (7) lists the Chandra ACIS Observation ID. 
Column (8) lists the observed X-ray luminosity in the 2-10 keV band, estimated as described in Section 2.1. 



WASPP pipeline identifies all stars in the USNO-A2.0 
Catalog that might be on the image, and searches for 
correspondence between the catalog and observed bright 
sources. The results are used to improve the astrome- 
try of the association data products. In principle the fi- 
nal astrometric accuracy could be as good as that of the 
USNO-A2.0 Catalog, which has a nominal RMS accuracy 
of 0.25". In practice, however, there are several caveats 
to this. First, the USNO-A2.0 Catalog has epochs in the 



1950s, and acciimulated proper motions bring the cur- 
rent RMS accuracy into the range 0.3"-0.4" (Schade et 
al. 2002). Second, not all fields have (identifiable) USNO- 
A2.0 stars, and for these fields the association data prod- 
ucts have the same 1.6" RMS astrometric accuracy as 
regular pipeline-calibrated WFPC2 data. And third, the 
nearby galaxies of interest for the present study generally 
fill the field of view and have complicated morphology. 
In this case the automated source detection software in 
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WASPP may identify the wrong source with a USNO- 
A2.0 star, and thereby significantly degrade the astroni- 
etry of the image. In our sample we found several cases 
for which this appeared to have happened. Many galaxy 
nuclei are themselves contained in the USNO-A2.0 Cat- 
alog, and due to their complicated nuclear morphologies 
these are particularly tricky to use in any astrometric 
registration. 

Because of these partial shortcomings of the WASPP 
astrometric registration algorithm, we did not use its 
results. Instead, we visually identified stars from the 
USNO-Bl.O Catalog (Monet et al. 2003) in the associ- 
ation data products, and used them to improve the as- 
trometry. While USNO-A was a two-color, one-epoch 
catalog, USNO-B instead is a three-color, two-epoch cat- 
alog. The main advantage in the present context is that 
proper motions are included in the catalog, and can hence 
be corrected for. The nominal RMS accuracy of the cat- 
alog is 0.2". Note that mosaics produced using the STS- 
DAS task wmosaic were used, which removes known rela- 
tive rotations and geometric distortions of the individual 
WPFC2 chips. The relative astrometry between WFPC2 
positions should then be good to better than 0.1" ^. The 
advantage of our visual identification of the stars is that 
it avoids the potential errors that an automated algo- 
rithm could make. This approach allowed us to correct 
the astrometry of 34 of our 76 images. For the remain- 
ing 30 images there were no usable USNO-Bl.O stars in 
the field (and in general, the WASPP pipeline had not 
identified any USNO-A2.0 stars either). In these cases 
we left the astrometry as contained in the header of the 
original (pre-association) lead image which came out of 
the STScI on-the-fly-recahbration (OTFR) pipehne. 

For the purpose of the present study we are not really 
interested in the absolute astrometry of either the Chan- 
dra or the WFPC2 data itself, but rather in the relative 
astrometric accuracy between the two. It is possible to 
increase this accuracy in those cases where the WFPC2 
field contains a known X-ray source that also emits at 
optical wavelengths (other than the ULX itself). This is 
the case if the galaxy nucleus is a point-like AGN, and if 
the nucleus is not obscured by dust in the WFPC2 image. 
We found this to be the case in 27 of the images, corre- 
sponding to half of all the galaxies in the sample. This 
high fraction reflects the fact that the Chandra and/or 
HST observers tend to bias their observations towards 
active galaxies. 

With the aforementioned approaches we obtained three 
types of WFPC2 header coordinates: Method 1 coor- 
dinates (available for 27 images) were shifted to align 
the galaxy nucleus with the AGN coordinates in Chan- 
dra data; Method 2 coordinates (available for 46 im- 
ages) were shifted to align observed USNO-Bl.O stars 
with their catalog coordinates; and Method 3 coordi- 
nates (available for all 76 images) contain the original 
astrometry from the telescope. These methods are listed 
in order of decreasing accuracy. For our final analysis, 
we adopt for each image the coordinates from the lowest 
method number available, as listed in Table 2. Out of 76 
images total, 27 images have coordinates from method 1, 
34 images have coordinates from method 2, and 15 im- 
ages have coordinates from method 3. Table 2 also lists 

3 

see http://www.stsci.edu/instruments/wfpc2/Wfpc2_faq/wfpc2_ast_faq.html 



for each image the final pixel position of each ULX. To- 
gether with the RA and DEC values of the ULXs listed 
in Table 1 this uniquely defines the astrometry of each 
image. Note that our overall approach focuses on relative 
astrometry and not absolute astrometry. In particular, 
when we register WFPC2 data to the Chandra frame in 
method 1 we improve the relative astrometry, but not 
necessarily the absolute astrometry. 

The relative accuracy of the coordinates from the dif- 
ferent methods can be analyzed with some simple statis- 
tics. Let dij be the RMS (two-dimensional) distance be- 
tween the coordinates from methods i and j, for those 
images that have coordinates available from both meth- 
ods. We find for our sample: di2 = 0.82" ± 0.11", 
di3 = 1.67" ± 0.16" and ^23 = 1-40" ± 0.13". We can 
safely assume that method 1 yields much better relative 
Chandra/HST alignment than method 2. This implies 
that di2 is a measure of the relative RMS Chandra/HST 
alignment for those images for which the best coordinates 
are available from method 2. Similarly, dis is a measure 
of the relative RMS Chandra/HST alignment for those 
images for which the the best coordinates are available 
from method 3. In the following we use di2 and di3 to set 
the size of the error circles in our search for ULX coun- 
terparts. For method 1 the relative RMS Chandra/HST 
alignment is more diSicult to assess, because it is domi- 
nated by the accuracy and stability of the geometric dis- 
tortion solutions for the two instruments. For WFPC2 
the relative positions of the chips have shifted by ~ 0.2" 
since launch (Anderson & King 2003). Furthermore, this 
alignment depends on the accuracy to which a centroid of 
a source can be determined in both the X-ray and the op- 
tical data. The complexity of the morphology can affect 
this determination. In the following we conservatively, 
and somewhat arbitrarily, assume a RMS Chandra/HST 
alignment accuracy of 0.3" for method 1. 

We note that our results are more-or-less consistent 
with expectation. For example, if the absolute RMS ac- 
curacy of USNO-Bl.O registered coordinates is 0.4", the 
absolute RMS accuracy of Chandra pipeline coordinates 
is 0.7", and the absolute accuracy of WFPC2 pipeline 
coordinates is 1.4", then we would have expected that 
di2 = 0.81", di3 = 1.57" and ^2.3 = 1-46". This is in sta- 
tistically acceptable agreement with our sample statis- 
tics. Absolute accuracies of 0.4" for USNO-Bl.O regis- 
tered coordinates and 0.7" for Chandra pipeline coor- 
dinates are slightly larger than the nominally expected 
values. However, it is not unreasonable to expect that 
a variety of small systematic errors could increase the 
achieved accuracies so as to exceed the nominal values. 

4. IMAGES 

Figure 1 shows large scale views of each of the galaxies 
and ULXs in our sample. Thc;se images are intended to 
illustrate the overall morphologies of the sample galax- 
ies and the relative locations of the ULXs with respect 
to the galaxy nucleus. We therefore show only data 
for a single filter per galaxy. When images in multiple 
filters were available we preferentially chose to show a 
broad- or medium-band filter in the red or visual part of 
the spectrum, such as F814W, F791W, F606W, F555W, 
or F547M (see Biretta et al. (2001) for a discussion of 
the properties of WFPC2 filters). However, there are 
some exceptions to this. For NGC 1097, NGC 5033, and 
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NGC 7174 there is only data in the uhra-violet (F218W, 
F300W) part of the spectrum. For the Antennae, one of 
its ULXs (#25) is included only in a U-band (F336W) 
image. For NGC 4594 only a narrow-band filter image is 
available (F658N). For Circinus, one of its ULXs (#41) 
is included only in a narrow-band image (F656N). 

The greyscale panels in Figure 2 show 10 x 10 arc- 
second regions centered around the individual ULXs 
in the sample. A circle around each ULX indicates 
the la accuracy of the astrometric registration, i.e., 
0.30" for images aligned by method 1, 0.82" for im- 
ages aligned by method 2, or 1.67" for images aligned 
by method 3. These images show the general galaxy 
morphology around each ULX. They form the basis for 
the search for optical counterparts described in Section 5. 
Figure 2 shows each ULX, but not each dataset for each 
ULX. For conciseness we generally only show the high- 
est quality image for each ULX. In most cases the mor- 
phology is visually very similar in the different filters for 
which data are available. However, in some cases the im- 
ages in different filters do show different morphologies, 
and for these ULXs we show images of the same region 
in multiple filters in Figure 2. 

The CADC WASPP pipehne performs a full reduction 
and image combination for the datasets in each associ- 
ation. We therefore did not attempt any additional re- 
duction or calibration of the images. However, we did 
discover occasional peculiarities in the combined images 
from the CADC WASPP pipeline. In particular, some of 
the images show evidence of improperly removed cosmic 
rays. Most of these are cases in which only two expo- 
siires arc available, and for which the data may simply 
having been inadequate to fully remove all cosmic rays 
and their coincidences. Examples of galaxies for which 
some datasets appear to have residual cosmic rays are 
those of, NGC 4594, NGC 5018, and NGC 5033. All 
counterparts were visually inspected to screen question- 
able cases. 

Another peculiarity in some of the images is a sprin- 
kling of negative pixels that seems inconsistent with the 
regular tail of the noise distribution. We learned (CADC, 

priv. comm.) that these are generally the result of com- 
bining two or more images with significantly different ex- 
posure times. We made sure that these pixels generally 
do not affect our photometry of potential counterparts 
(unless in a few rare cases, which we explicitly identify 
as such below). Figure 2 also shows some peculiarities 
that are intrinsic to the data and the WFPC2 instru- 
ment. Most notably, the diagonal bright stripes visible 
in the images for ULXs 3,4, and 31 are chip boundaries 
in the mosaiced data. 

5. PHOTOMETRY 

We searched all the images for potential ULX optical 
counterparts. In doing so we restricted our attention to 
a search radius in the optical images that is twice the as- 
trometric registration error for each ULX. All the ULXs 
in our sample reside in nearby galaxies, and many of 
these galaxies are spiral galaxies. Such galaxies can have 
complicated morphologies, and we therefore commonly 
found diffuse structures of various kinds (e.g., parts of 
spiral arms, dust bands, etc.) within the search radius 
around a ULX. Such structures may be interesting if a 
connection to the ULX phenomenon could be demon- 



strated. However, we decided not to catalog such struc- 
tures in the present context. Instead, we decided to re- 
strict our attention to sources within the search radius 
that appear point-like, i.e., sources that are either consis- 
tent with the HST point spread function (PSF) or which 
are only slightly more extended. Such sources tend to 
be associated with star clusters, diffuse regions, or indi- 
vidual stars. All of these can provide essential insight 
into the ULX phenomenon if a direct association with a 
particular ULX could be demonstrated. 

Our method for the identification of point-like sources 
depended on the galaxy morphology around the ULX. 
Based on visual inspection we discriminated between sev- 
eral different classes of morphologies within the search 
radius. We classify as simple morphologies those cases in 
which there are one or more easily identifiable, point-like 
sources on a relatively smooth or constant background. 
For these morphologies we used the standard source de- 
tection routine daof ind in the IRAF software package to 
find point-like sources. We utilized 4cr as the detection 
threshold. We will refer to datasets in this category as 
case S (S for "simple" ) . Other ULXs arc found in regions 
of more complicated morphology. These show diffuse, ex- 
tended or dusty structures. In these cases we found that 
automated routines for source detection did not provide 
reliable results. Nonetheless, in some fraction of cases 
with these complicated morphologies it was possible to 
manually identify point-like sources. We will refer to 
these as case M (M for "manual"). We emphasize the 
point that the source detection in these cases is never- 
theless somewhat arbitrary, particularly with regard to 
whether faint sources are considered to be real or arti- 
facts. For the remainder of the complicated morpholo- 
gies it was not possible to place any reliable limits on the 
presence or absence of point-like sources. We do not re- 
port any results for these cases, which we will refer to as 
case C (C for "complex"). This primarily occurred in the 
cases of two galaxies, NGC 1068 and NGC 3034 (M82). 
In the case of NGC 1068, while it appears that in some 
fields manual source detection would be possible, taken 
as a whole the combination of large numbers of sources 
within the 2a error circle (3.3" in radius) and complex 
spatial morphology would complicate the interpretation. 
We therefore defer detailed analysis in this case for future 
work when better astrometry may be available (e.g., by 
tying the astrometry of WFPC2 images to ground-based 
images). We also classified a source as "C" due to it 
being on the edge of the chip (a possible counterpart to 
ULX 24 in the u3040205b observation). 

Another category is formed by those ULXs for which 
the optical morphology within the search radius is 
smooth and straightforward, but in which there is no sign 
of detectable point-like sources. We will refer to these as 
case N (N for "nothing" ) . A final category is formed by 
those ULXs for which the optical data was unreliable due 
to a high fraction of unrealistically negative pixels within 
the search radius. This problem of the data calibration 
was discussed in Section 4. It only affected a few datasets 
so significantly that we felt that a quantitative analysis 
would not be warranted. These datasets are referred to 
as case X (X for "excluded"). For each ULX and dataset 
we indicate in column (5) of Table 3 to what case it be- 
longs. Of the 145 unique dataset, ULX, and photometric 
case combinations, 20 (14%) represent case S, 35 (24%) 
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Fig. 1. — WFPC2 mosaics showing the large scale galaxy morphology and ULX position for all of the ULXs in our sample. Each 
ULX is shown as a cross, marked by its ID number in Table 1. Each image is 2.58 x 2.58 arcmin, and has North pointing towards 
the top of the page. The galaxy name, dataset, and filter are indicated at the top of each panel. For a few galaxies we show multiple 
images for different pointings, because not all ULXs fall on a single image. Only first page of figures are shown here, for full version see 
http: / /xassist. pha.jhu.edu/ptaJi/hst_ulx_paper. 
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#1, NGC 278. u66Cie601b. F450W HQ!} 27S, u6eaeS05b. F814W 
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Fig. 2. — WFPC2 mosaics showing the 10"xlO" regions around the ULX positions with well-defined measurable counterparts. North is 
pointing up. The picture on the left shows the region of the galaxy with a one sigma error circle drawn on top. Tiic picture on the right 
shows both the one and two sigma error circles, and the positions and numbers of the optical counterparts ovcrlaycd and marked with 
crosses. For discussion of the error circle sizes, sec Section 3. Details about the photometry for the counterparts are shown in Table 3, and 
discussed in Section 5. Only first page of figures are shown here, for full version see http;//xassist. pha.jhu.edu/ptak/hst_ulx_paper 
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represent case M, 37 (26%) represent case C, 44 (30%) 
represent case N and 9 (6%) represent case X. 

For 9 ULXs we have one or more datasets of case S, 
and for 20 ULXs we have one or more datasets of case 
M. These are the cases for which one or more point-Hke 
sources were either automatically or visually identified 
within the search region around the ULX. In total, we 
found potential counterparts for 28 of the 44 ULXs in 
our sample ( there is overlap between S and M cases for 
ULX 40). We assigned each of the sources thus identified 
a unique ID number. The sources are labeled and shown 
schematically in separate panels of Figure 2. Informa- 
tion on the sources is provided in Table 3. Columns (7) 
and (8) give the offset of each source w.r.t. ULX position 
(in RA and Dec). Column (9) gives the correspond- 
ing total offset. For each source we performed aperture 
photometry using the routine phot in the IRAF software 
package. The photometry on the source was performed 
with an aperture of radius 0.1" (this is 1 pixel on the 
Wide Field Camera CCDs and approximately 2 pixels 
on the Planetary Camera CCD). An appropriate back- 
groTind value was subtracted from the results. For case S, 
the local background was automatically measured in an 
annulus surrounding the point source. For case M, the 
background was kept fixed at the value determined man- 
ually from a representative (but carefully chosen) region 
near the ULX. 

The measured count rates in apertures were converted 
to magnitudes using the zeropoints for each filter given 
in the WFPC2 Data Handbook (Mobasher et al. 2002). 
The resulting magnitudes are in the so-called VEGAMAG 
system. This is a synthetic system in which the mag- 
nitude of the star Vega is defined to be exactly zero in 
all filters. In the traditional Johnson-Cousins system the 
magnitude of Vega is approximately zero in the UBVRI 
bands, but not exactly (Bessel, Castelli, & Plez 1998). In 
general, magnitudes on the VEGAMAG system are close to 
those on the Johnson-Cousins system, provided the filters 
have similar central wavelengths. However, to actually 
transform the magnitudes onto the Johnson- Cousins sys- 
tem it is necessary to know the spectral shape, or at least 
a color, of each of the identified sources. Since this is not 
generally known, or known well, we have not attempted 
such transformations. 

Various corrections are necessary to the inferred mag- 
nitudes. First, an aperture correction is necessary to 
estimate the true magnitudes from those that are mea- 
sured in a small aperture. We used the TinyTim soft- 
ware (Krist & Hook 2001) to generate synthetic PSFs 
for each of the filters. From aperture photometry on 
these PSFs we calculated the aperture corrections be- 
tween the apertures that we used for the data analysis 
and a 0.5" radius aperture. These corrections were ap- 
plied to the inferred magnitudes. An additional correc- 
tion of 0.1 mag was applied to transform the results to 
a nominal infinite aperture (Mobasher et al. 2002). A 
second correction that must be applied is for geometric 
distortion, which results in a change of the surface area 
covered by each pixel. This pixel area variation is im- 
printed on point source counts during the fiat-fielding 
step (which is defined to make a source of uniform sur- 
face brightness look "fiat" on the detector). We corrected 
for this using the known pixel area variations (Mobasher 
et al. 2002). The correction amounts to at most 0.05 



mag in the corners of each chip. The third correction 
that must be applied is for imperfect charge-transfer ef- 
ficiency (CTE), which causes charge to be lost from the 
aperture during read-out. We corrected all photometric 
results (detections only) using the CTE correction al- 
gorithm of Dolphin (2000), using the equations posted 
at http://purcell.as.arizona.edu/wfpc2_calib/ (Decem- 
ber 2004 version). The CTE corrections increase with 
decreasing background. 

The final magnitudes obtained after application of all 
corrections are listed in column (10) of Table 3. We cal- 
culated the total statistical error as [eg + 2A{ag)^Y/'^ , 
where c is the net number of counts in the aperture, A 
is the area of the aperture in pixels, a is the noise level 
per pixel determined from the background region, and 
g is the gain (either 7 or 14 electron/data unit). This 
error is more conservative than the error computed by 
the phot routine and was on average 0.02/0.05 magni- 
tudes larger for sources brighter /fainter than 24. It is 
likely that these are underestimates of the true errors, 
especially in cases with complicated morphologies. Im- 
plied absolute magnitudes were calculated based on the 
galaxy distances quoted in Table 1 and are listed in col- 
umn (12). Signal-to-noise ratio (S/N) estimates were 
calculated for each counterpart in each image and are 
listed in column (13), using the formula S/N = ^^o.a ■ 
Sources for case S mostly had S/N > 4, due to the detec- 
tion threshold used in the daof ind routine. For case M 
sources we found that our manual source identification 
was also able to detect sources down to S/N « 4, al- 
though most of the sources thus identified have much 
higher S/N. 

For each ULX and each dataset (except those of case X 
for which the data were suspect and case C where the 
local background level is ambiguous) we estimated the 
limiting magnitude down to which a source might plau- 
sibly have been detected. For this we calculate the noise 
in a background aperture of the same size as was used 
for source photometry. We multiply the result by 4 to 
get the counts needed for a 4cr detection. We transform 
this result to magnitudes using the exposure time and 
filter zeropoints and apply the appropriate aperture cor- 
rection. The results are listed in column (6) of Table 3. 
For comparison, we also calculated the limiting magni- 
tudes using the Exposure Time Calculators (ETCs) for 
WFPC2 maintained by STScI. These more idealized cal- 
culations generally give somewhat fainter limits than the 
ones we calculated. This is as expected, given that the 
ETCs do not take into account that we are searching for 
sources in nearby galax;ies. The galaxy light adds both 
noise, fluctuations and complexity. 

In the final column we list source radius estimates. We 
found that kernel-fitting routines such as imexamine did 
not converge consistently, particularly for "M" sources. 
We therefore estimated the source sizes visually from the 
mosaic images (with a uniform pixel size of 0.1"). While 
siibjective, the intent here is to simply give an estimate 
for the spatial extent of the sources and not a detailed 
spatial analysis. We plot in Figure 3 a histogram of 
the estimated source sizes. The majority of sizes were 
in the range 0.1-0.4", with a few outliers having larger 
sizes. Excluding these outliers, the mean and standard 
deviation of the source radii are 0.22" and 0.05" , respec- 
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Fig. 3. — The size distribution of the ULX counterparts. Tiie 
sizes plotted are radii estimated visually from the mosaic images. 
These source sizes are probably representative of point sources, 
where the mean radius corresponds to ~ 85% of the encircled en- 
ergy fraction (see text for details). 

tively. For comparison, from the WFPC2 Instrument 
Handbook, this mean radius corresponds to an encircled 

energy fraction of 85% at 4000A. In cases where the 
use of imexainine resulted in a reasonable full-width half- 
maximum (FWHM) estimate, the mean value was 0.15" 
whereas the mean visual source size for these sources 
was 0.19". Therefore most sources are likely to be point 
sources or at most only marginally resolved, with our vi- 
sual size estimate being ~ 30% larger than the FWHM. 
We also list the linear extent corresponding to the source 
size assuming the distances listed in Table 1. 

6. DISCUSSION 

We have presented a systematic analysis of archival 
WFPC2 observations for a sample of ULXs. The main 
motivation was to determine the characteristics of possi- 
ble optical counterparts to ULXs. Minimally a key goal 
was to determine what is the probabihty of detecting 
a counterpart. Other goals would be to establish the 
typical optical luminosity of cotmterparts and how the 
implied X-ray to optical (or UV) luminosity ratio com- 
pares to known source types. Finally, in cases in which 
observations are available in multiple wide-band filters, 
the optical colors might also be a useful discriminator of 
sotirce type. Below we address the color information con- 
tent of the sample prior to computing optical luminosities 
since any knowledge of the spectral types of the counter- 
parts would benefit the computation. Of course, when 
multiple counterparts are present, obviously at most only 
one of the counterparts within the error circle can be the 
true counterpart of the ULX. Also, when only a single, 
faint counterpart is observed, particularly in late-type 
galaxies, the association may simply be due to chance 
coincidence. Therefore this discussion should be taken 
as an assessment of what properties any potential coun- 
terparts would have on average. 

6.1. ULX WFPC2 Counterpart Frequency 

As is evident from Table 3 and Figure 2, variably zero, 
one or many counterparts are detected within the (com- 



bined X-ray/HST) error circles. Since the stellar popu- 
lations in elliptical galaxies are old (and hence faint), a 
simple assumption would then be that finding no coun- 
terpart would be more common in early-type galaxies. 
Also, since there is an association between X-ray bi- 
naries and globular clusters in elliptical galaxies (see, 
e.g., Maccarone, Kundu, & Zcpf 2003), we might ex- 
pect to find counterparts to extra-nuclear X-ray sources 
in elliptical galaxies that are likely to be globular clus- 
ters. To test these assumptions, we plot in Figure 4 his- 
tograms showing the Hubble type distribution for the 
Chandra/WFPC2 galaxy sample (with 25 galaxies to- 
tal), the subset in which at most one counterpart was 
detected (15 galaxies) and the subset in which more than 
one counterpart was detected and/or an observation was 
marked as complex (10 galaxies). As expected, there 
is a tendency for single or no counterparts to be found 
in early- type galaxies and multiple counterparts and/or 
complex cases to occur in late-type galaxies. 

In Figure 4 we also compare the Hubble types of the 
galaxies in our Chandra/WFPC2 sample to the type dis- 
tribution of nearby (cz < 5000 km s~^, or D < 66 Mpc 
for Ho — 75 km s~^ Mpc^^ and go = 0.5) galaxies in 
RC3, in order to assess if there is any bias in our galaxy 
sample. While the statistics are limited, there is a slight 
excess of ellipticals and irregular/peculiars, and a more 
pronounced excess of Sb/c galaxies in our sample rela- 
tive to nearby galaxies. Note that the prominent peak in 
the Chandra/WFPC2 galaxy type distribution at Sb/c 
is mainly due to a large number of Sb (type = 3) galaxies 
(see Table 1). Therefore, in addition to the excess of Sb 
galaxies, the main bias in the sample is a lack of early- 
type (i.e., SO, Sa) and late type (i.e, Sd) spiral galaxies. 

Finally, we computed the number of background 
sources expected for the galaxies with Hubble type <3 
and Hubble type > 3. This was done by taking the larger 
of the limiting flux in each Chandra field and the flux as- 
sociated with a Lx = 10"^^ ergs s^^ source (i.e., see Ptak 
& Colbert 2004), in both cases computed in the 2-10 
keV band based on the full-band count rates. The num- 
ber of background sources expected in each field was then 
taken from the 2-10 keV logN-logS plot given in Bauer et 
al. (2004) derived from the Chandra Deep Fields. This 
resulted in an estimate of 5.5 background sources in the 
early-type galaxies (35% of the total) and 3.8 background 
sources in the late-type galaxies (14% of the total). 

6.2. Optical Colors 

We turn now to cases in which there are observations of 
a ULX in more than one wide-band filter and at least one 
detection (and hence a color or color limit can be com- 
puted) . Here we restrict our sample to observations with 
the filters F336W, F450W, F555W, F606W and F814W 
since they map reasonably well to standard UBVI mag- 
nitudes. With both F555W and F606W approximating 
the V band, hereafter we refer to them as V555 and Veoe- 
Only a few ULXs did not have a (CADC) observation 
in one of these filters, and in those cases the observa- 
tions were entirely in either narrow-band (e.g., ULX 41 
in Circinus) or UV (e.g., ULXs 43 and 44) filters. The 
colors and color limits are listed in Table 4. The colors 
B-V and V-I are plotted (as a function of Hubble type 
of the host galaxy) in Figure 5. Note that 6/9 of the 
Veoe-I and 5/8 of the B-Veoe data points are due to ULX 
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Fig. 4. — The Hubble type distributions of galaxies in the Chandra/WFPC2 sample, (left) The full (blue, solid line) sample distribution 
along with the subset in which at most one counterpart was detected (black, dashed line) and the subset in which more than one counterpart 
was detected and/or an observation was marked as complex (green, dotted line), (right) The fractional type distribution of galaxies in the 
full Chandra/WFPC2 sample (blue, dashed line) along with the type distribution of nearby {cz < 5000 km s~^) RC3 galaxies (black line). 



1 in NGC 278 (Hubble type 3 = Sb). The color F555W 
- F814W differs from Johnson V-I negligibly for V-I in 
the range -1 to 3 (see Holzman et al. 1995 '^). The sense 
of the correction between Veoe and V555 is such that it 
is negligible for sources with V-I < 0.5. For red sources, 
^555 - VeoG 0.3 - 0.6 (for stellar types G and MOIII), 
i.e., at worst we are underestimating V-I where V=V606- 

6.3. Optical Luminosities 

Excluding the counterparts marked as complex, there 
are 67 unique counterparts in which one of UBVI was 
used. For observations in these filters we first com- 
pute the absolute magnitude and the luminosity in solar 
units, e.g., logiy = — 0.4(Afy — Mq v), where we use 
Mq^u = 5.51, Mq,s = 5.41, MQy = 4.79, Mqj = 4.03 
for solar absolute magnitudes (Zombeck 1990). The re- 
sults are listed in Table 5. Since ^ 50% of the counter- 
parts include either a V555 or Vgoe measurements, and 
V is a convenient magnitude for comparison with pub- 
lished photometry and models, we proceed by comput- 
ing an effective V-band luminosity. For cases with color 
information available, we classified the counterpart as 
red if V-I > 0.5, B-V > 0.5, B-I > 1.0, or U-V > 0.7 
and blue otherwise. If no color information is available, 
then we assumed sources in galaxies with Hubble type 
< 3 to be red and blue otherwise. For red (blue) coun- 
terparts, we computed assuming a G5 (AO) stel- 
lar spectrum. We then computed Mv,eff, the effective 
absolute V magnitude, and Ly = vL^, (we assumed 
an effective wavelength of 5500A and zero-point flux 
density of 3.67 x 10^^*^ ergs cm^^ s^^ Hz^^; Zombeck 
1990), also listed in Table 5. Ly is listed in units of 
Lq^ = 2.9 X 10'^^ ergs s~^, and these values are in good 
agreement with those computed from the absolute mag- 
nitudes in individual filters. 

In Figure 6 we show the distributions of Veff and 
logiy. While the statistics are limited, the Veff dis- 



tribution is reminiscent of a power-law distribution with 
a cut-off at y ~ 24, which is consistent with the average 
hmiting magnitude of the fields discussed here (24.9). 
The logLy distribution shows that we are observing 
sources at least as luminous lO^Lo^y, with a mean value 
of ~ IQ^L^y. We also plot the luminosity distributions 
only including sources brighter that V=24 and V=23.5. 
The large drop-offs in number counts below ~ IO^Lv^q 
suggests that we are only approaching completeness for 
counterparts with L > IQ^Ly^Q- We also plot the effec- 
tive V band luminosity distribution for upper-limits (also 
listed in Table 5). For the majority of cases the limiting 
luminosity was 10^ or higher, suggesting that the lack 
of a counterpart was due at least in part to the lack of 
sufficient depth in those fields to detect counterparts of 
comparable luminosity to the counterparts in the entire 
sample. 

An important issue is to what extent there is a trend 
in counterpart luminosity with distance since at larger 
distances only brighter sources would be detectable (i.e., 
the Malmquist bias), and source confusion may also be- 
come problematic at larger distances. In Figure 7 we plot 
the the optical luminosities of the counterparts as a func- 
tion of distance, and there is in fact a strong correlation 
between logLy and distance, although only at distances 
< 15 Mpc. We also plot the luminosity expected for a 
source at a constant fiux level (with arbitrary normal- 
ization) to show the sense of the Malmquist bias. The 
observed optical luminosities increase more rapidly than 
this trend, although of course in practice the observa- 
tions were performed with a range of filters and exposure 
times resulting in a large range in limiting flux. Source 
confusion may also be contributing, however the lack of 
a strong trend at larger distances suggests the neither 
source confusion nor Malmquist bias is a factor at dis- 
tances > 20 Mpc. Since our discussion is largely based 
on optical luminosities spanning ^ 2 dex, the presence 
of some source confusion is not going to strongly impact 
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Fig. 5. — V-I (left) and B-V (right) colors for counterparts to ULX sources. F891W and F450W have been used as proxies for I and B, 
while F555W and F606W have been used as a proxy for V (plotted with black and red points, respectively). Note that a large fraction 
of the data points in these plots are due to counterparts for ULX 1 (see Table 3). The shaded regions show the l-cr range in global color 
determined for spiral galaxies in de Jong (1996). For reference, the WFPC2 Exposure-Time Calculator Tool was used to compute colors 
using the Sbc galaxy spectral model (T=4), which resulted in B-V = 0.4 (0.6) and V-I=1.0 (0.8) for V=F555W (F606W), consistent with 
the de Jong values. 
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Fig. 6. — Histograms of the computed effective V magnitudes (Ve//i left) and the corresponding luminosities in solar units (right). The 
luminosity plot also includes the cases of only including sources brighter than Vgff = 23.5 (dashed line) and Vgff = 24 (dot-dashed line), 
as well as upper-limits (dotted line). This shows that this sample is only approaching completeness for Lv/Lq.v > 10®. 



our conclusions. 



6.4.1. Isolated X-ray Binary 



6.4. Source Assessment 

Unless an ULX is in reality an interloper (i.e., a fore- 
ground star or background AGN), the most likely phys- 
ical model is that they are some sort of X-ray binary 
(XRB; e.g., see Colbert et al. 2004). In this case, the 
most optimistic scenario is that a given optical counter- 
part is the companion star. However, the ULX coun- 
terpart may also represent a cluster of stars or simply a 
clump of gas or stars (particularly in an arm of a spiral 
galaxy). We address these in turn below. 



The absolute magnitudes in Table 5 are in the range -5 
to -12, corresponding to logLy/LQ^ ~ 4 — 6. Blue gi- 
ant stars have log Ly / Lq^ ^4 — 5 and supergiants can 
have logLy/ ~ 5 — 6. In the supergiant case the 
colors are not restrictive since blue, yellow or red super- 
giants could be present. However in all of these cases 
the stars would be very young (< 10^ years; Bertelli 
et al. 1994) and massive (> IOM0) and therefore the 
XRB would be a high-mass XRB. High-mass XRB typi- 
cally have relatively low X-ray/optical flux ratios, in the 
range of -2 < logL2-ii kcv/isooo-Tooo < 1 (Bradt & 
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Fic;. 7. — Effective V-band luminosity plotted as a function of 
distance to the host galaxies. The line plotted shows the trend 
expected for a source with constant flux (i.e., Ly oc D^, where D 

is the distance). 

McClintock 1983). In Figure 8 we show the distribution 
of log Lx/Lv.eff values from our sample along with the 
observed values of log_L2-ii kcv/-^30oo-7ooo from Bradt 
& McClintock (1983) for low and high-mass XRB. In gen- 
eral, \ogLx/Lv,eff ~ logL2-ii kev/^aooo-Tooo- By def- 
inition ULXs are orders of magnitude brighter than typi- 
cal high-mass XRB however the stellar companion would 
not necessarily be more luminous than those found in 
Galactic high-mass XRB. Accordingly the \ogLx / Lv,ef f 
values for our counterparts would be expected to be 1-2 
dex higher than those found high- mass XRB, and Figure 
8 is consistent with this. Therefore the high-mass XRB 
scenario is not ruled out by our data, although of course 
such young stars would not be expected in early-type 
galaxies and are therefore unlikely in galaxies with Hub- 
ble types less than ~ or for ULXs that are clearly bulge 
or halo sources, such as ULX 32 in NGC 4565. Con- 
versely a high-mass XRB association with ULXs would 
be consistent with the ULX/star formation connection 
that has been suggested (e.g., Swartz et al. 2004). 

Another possibility is that the ULX is an X-ray binary 
with an accretion disk, and that the accretion disk is 
contributing significantly to the optical-UV part of the 
spectrum. Copperwheat et al. (2005) discuss this possi- 
bility, taking into account the impact of X-ray irradiation 
of the accretion disk and companion star. They find that 
the expected V-band absolute magnitudes {My) would 
be in the range of -4 to -9, with the companion star 
often dominating over the accretion disk at low black 
hole masses (10 < Mbh < lOOAfg) luilcss irradiation 
is significant. This implies optical counterparts with 
3.5 < logLv < 5.5, consistent with many of our counter- 
parts (see Figure 6). Copperwheat et al. also generally 
predict blue colors for the disk + star system. While 
those results were calculated for the case of an isotropic 
X-ray luminosity of lO'*" ergs s~^, we find similar results 
for higher luminosities (and hence higher accretion rates 
for a given black hole mass). Note that the assumption 
of the X-ray luminosity being 10% of the bolometric lu- 
minosity gives log Lx/Lv ~ 1-0 — 1.6, which is often 
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Fig. 8. — The distribution of log Lx / Ly^f f from our sam- 
ple. Also shown arc the distributions in logL2— 11 kev/^3000— 7000 
for low-mass (dotted line) and high-mass (dot-dashed line) X- 
ray binaries in Bradt & McClintock (1983). The correction from 
log Lx/Ly^ff to L2-11 kcV log ^3000-7000 should be small since 
the peak wavelength of the V band, ~ 5500A, is near the midpoint 
of 3000-7000A. 

consistent with our observed values. Finally, similar re- 
sults to those discussed above would be expected in the 
case of Bondi accretion onto an intermediate-mass black 
hole in a molecular cloud (Krolik 2004), with the optical 
emission originating several hundred gravitational radii 
from the black hole as in the thin-disk case (J. Krolik, 
priv. comm.). 

6.4.2. Structure in Spiral Disks 

Inspection of Figure 2 suggests that in several cases 
we are observing fine structure in the spiral arm of the 
galaxy. In this case the counterparts would likely be local 
maxima in the surface brightness and/or holes in the 
extinction. The mean source extent of 0.2" corresponds 
to a linear size of ~ 20 pc for the mean distance of 20 
Mpc for the galaxies in this sample (see Table 3). In 
Figure 5 we also plot the l-cr range in B-V and V-I color 
for spiral galaxies given in de Jong (1996). The galaxies 
NGC 278 and NGC 891 are the only spiral galaxies (both 
Sb galaxies) with color information for counterparts, and 
again most of the data points are due to ULX 1 in NGC 
278. While the data are (very) limited, the colors are 
bluer than expected for the integrated emission of an Sb 
galaxy, possibly implying regions that are younger than 
Sb disks are on average. This also would be crudely 
consistent with the ULX/star formation connection that 
has been suggested. 

6.4.3. Clusters 

We next consider the case of the counterparts being 
stellar clusters. In the case of early-type galaxies this 
would of course correspond to globular clusters (GCs). 
Although the numbers are limited, for earlier type galax- 
ies, B-V and V-I tend to be red (> 0.5), possibly indicat- 
ing an older stellar population (e.g., age > 10^~^ yrs; see 
Leitherer et al. 1999; Bruzual & Chariot 2003). For com- 
parison, the mean B-V and V-I colors for GCs in M31 
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(the Milky Way) arc 0.72 (0.71) and 0.96 (0.94) (Barmby 
et al. 2000), and the globular clusters in NGC 1399 asso- 
ciated with ULXs 8 and 9 (see Angelini et al. 2000) have 
B-V = 0.5, 0.9 and V-I=l.l, 1.0. Chandar, Whitmore, 
& Lee (2004) present an HST survey of GCs in nearby 
spiral galaxies and find that the GC colors similarly lie in 
the range 0.5 < B-V < 1.2 and 0.7 < V-I < 1.5, and that 
the GC luminosities peak at My ~ —6 to —8, consistent 
with most counterparts in our sample. The ULX coun- 
terparts in NGC 3311 (a lenticular) and ULX 11 in NGC 
1700 (an elliptical) have red colors (V-I > 1-2), and ULX 
40 in Ccn-A has a counterpart with V-I ^ 0.9. These arc 
likely to be GCs also. The ULX counterpart in NGC 
4565 has also been identified as a globular cluster (albeit 
with a blue color) by Wu et al. (2002). 

Conversely, V-I and B-V < 0.5 would be indicative of a 
stellar population that is very yoimg (< 10^ yrs) and/or 
has low metallicity (Leitherer et al. 1999; see also the 
discussion in Soria et al. concerning the application of 
stellar evolutionary models to counterparts in the vicin- 
ity of a ULX in NGC 4559). V-I and/or B-I colors in 
the range 0.0-0.5 are observed for several coimterparts 
in late- type galaxies (ULX 2 in NGC 891 and several of 
the counterparts for ULX 1 in NGC 278). Only U-V is 
available (i.e., with our CADC Association data selection 
criterion) for the ULXs in the Antennae. These counter- 
parts are all blue, with U-V < 0.5, and the coimterparts 
in the Antennae that appear more red are probably in 
regions of higher extinction. Also note that U-F336W 
varies from to 0.8 as U-B varies from to -1.5 (Holtz- 
man et al. 1995), and therefore the very blue U-V colors 
in Table 4 are over-estimated (but the true U-V are nev- 
ertheless blue). Interestingly there are also some coun- 
terparts in elhptical galaxies with blue colors (e.g., ULX 
37 in NGC 5018 with U-V = -1.5). If associated with 
the galaxies these are likely to be young stellar clusters 
captured or formed after a recent encounter or merger. 
Finally, there are counterparts in late-type galaxies with 
red colors (specifically several of the counterparts to ULX 
1 in NGC 278 and ULXs 15-16 in NGC 3256) however 
this may be due to extinction in these galaxies. Alterna- 
tively the red colors may be due in part to the onset of 
a red supergiant phase, although this would require the 
rather specific circumstances of a ^ 10^ year old clus- 
ter resulting from an instantaneous star formation burst 
(Leitherer et al. 1999). 

Even if the optical counterpart is a cluster, the X-ray 
emission is still likely to be due to an accreting X-ray bi- 
nary. In this case the X-ray/optical flux ratio would be 
"diluted" by the optical emission of the stars from the 
cluster, i.e., the observed X-ray/optical flux ratio would 
be lower than that of an isolated X-ray binary, by a fac- 
tor of Lv.c/ Ly.b — 1, where Ly.c is the cluster V-band 
luminosity and Ly.t is the X-ray binary V-band luminos- 
ity. As discussed above, some of the counterparts in our 
sample have log Ly > 5.5, in excess of what might be ex- 
pected from an isolated ULX (Copperwheat et al. 2005) 
by at least 1-2 dex. In the case of high-mass X-ray bina- 
ries (HMXB), the distribution of the expected Lx/Ly 
ratio largely overlaps the observed distribution of our 
counterparts, although the peak of the HMXB Lx/Ly 
distribution is ~ 2 dex lower than the peak of the coun- 
terpart Lx/Ly. Therefore the counterpart Lx/Ly dis- 
tribution is consistent with HMXB distribution assuming 



on the order of 1% of the cluster optical emission is due to 
the HMXB. Conversely, in old clusters (> 10^ years), we 
would expect any X-ray binary to be a low-mass X-ray 
binary (LMXB). Here larger amounts of dilution would 
be consistent with the observed Lx/Ly^^ff values since 
the typical LMXB Lx/Ly^eff values are 2 dex or more 
larger than observed in our counterparts. In summary, 
our observed Ly and X-ray/optical flux ratios are gener- 
ally consistent with a LMXB in old clusters, but probably 
only consistent with a HMXB in young clusters if the op- 
tical light of the X-ray binary is a larger fraction of the 
optical light of the entire cluster than is the case with a 
LMXB in an old cluster. 

Finally, we estimate the black hole mass that would 
be expected if a black hole were present with ~ 0.1% 
of the cluster mass. This is roughly consistent with 
the well-known bulge mass / black-hole mass trend, al- 
though the slope is highly uncertain at the low mass end. 
Laor (2001) finds that Mbh iM^uige tends to be lower for 
lower-mass systems (i.e., MBH/Mi,uige ~ 0.5% in ellip- 
ticals versus MBH/Mhuige ~ 0.05% in late-type spirals). 
However recent modeling has shown that clusters may 
form a black hole with mass ~ 0.1% of the cluster mass 
in a relatively short time (< 10 Myr; Portcgics Zwart & 
McMUlan 2002; Geurkan, Freitag, & Rasio 2004). The 
results are shown in Figure 9 and there is no obvious 
trend between Ly (and hence estimated Mbh) and Lx- 
Note that we plotted black hole mass estimates based on 
a mass-to-light (M/L) ratio of 1, i.e., Mbh = ^ogLy — 3., 
which is only appropriate for a cluster of age several bil- 
lion years. Prom population synthesis modeling the M/L 
ratio of a cluster is expected to range from 0.05 to 5 
for clusters of age lO"^ to 10^° years (McLaughlin & Van 
der Marel 2005). Therefore for young clusters the black 
hole masses would be ~ 1 dex lower than plotted here. 
We also plot the Eddington luminosity for a given black 
hole mass, LEdd = 1-3 x 10^^ ergs s~-^, which is the 
maximum expected X-ray luminosity for isotropic emis- 
sion. Most of the X-ray luminosities are well below the 
predicted Eddington luminosity given the estimated Ly, 
although there are some cases where the X-ray luminosity 
exceeds this value, suggesting either anisotropic emission 
or a larger black hole mass than ~ 0.1% of the cluster 
mass (again the M/Ly ratio may differ significantly from 
1 as assumed here). Also note that, particularly in the 
case of spiral galaxies, multiple counterparts are found 
for some ULXs, and the corresponding BH masses for 
each counterpart are plotted in Figure 9. Of course it is 
unlikely that all counterparts in those cases are bound 
systems with an intermediate-mass black hole. 

6.4.4. Interlopers 

We turn now to the case of either a foreground star 
or a background AGN. As stated above, we calculate 
that ~ 9.3 of the ULXs are likely to be background or 
forcgroimd X-ray sources (^ 21%). The typical values of 
X-ray to optical flux ratio expected in the case of an AGN 
would be — 1 < log Lx/Lopi ^ 1 (see, e.g., Norman et al. 
2004; Bauer et al. 2004) while in the case of foreground 
stars it would be < — 1 (c.f., Georgakakis et al. 2004). 
From Table 5, most of the counterparts are consistent 
with the background AGN possibility, particularly in the 
early- type galaxies, although a few have very low or very 
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Flc;. 9. — Estimated black hole mass assuming the optical coun- 
terparts are clusters with 0.1% of their cluster mass in a black hole 
(and a mass-to-light ratio of 1) plotted as a function of the ULX 
X-ray luminosity. More generally the M/L ratio of a cluster ranges 
from ~ 0.05 to ~ 5 for ages from 10^ to 10"' years (McLaughlin 
& Van der Marel 2005), and the black hole mass estimate plotted 
above would scale accordingly. The line shows the black hole mass 
expected if the ULX were emitting at the Eddington luminosity in 
X-rays. Note that in many cases there are multiple counterparts 
for a ULX. 

high X-ray to optical ratios, making a foreground star 
or an X-ray binary in one of the scenarios above more 
likely, respectively. 

7. CONCLUSIONS 

We have presented a comprehensive analysis of the 
available WFPC2 Association data with coverage of a 
sample of ULXs. The relative astrometry between the 
Chandra and HST data was corrected whenever possi- 
ble, in some cases requiring shifts of several arcseconds. 
The final error circles were 0.3-1.7" depending on the 
availability of sources useful for the registration. With 
the error circles there were varyingly no, one or many 
potential counterparts. The limiting luminosities are of- 
ten above the faintest observed luminosities suggesting 
that the lack of multiple counterparts in some cases is 
due to insufficient depth. We focused on the "wide" fil- 
ter observations corresponding to UBVI filters since the 
colors could be compared to known source properties. 
This analysis showed that there was a slight tendency 
for early-type galaxies to have counterparts that were 
red (i.e., B-V > 0.5, V-I > 0.5) and conversely late-type 
galaxies had blue counterparts. In general the source 
properties did not differ significantly from the properties 
of those expected from either structure in spiral disks 
(although in a few cases the colors were blue suggesting 
a young region) or clusters. The counterparts found in 
early-type galaxies are generally consistent with globu- 
lar clusters in both their colors and luminosities. The 
X-ray /optical flux ratios were generally consistent with 
high- mass X-ray binaries, background AGN or low- mass 
X-ray binaries in clusters (i.e., with artificially lowered 
flux ratios). Therefore for most of the sources, any one 
scenario cannot be ruled out although there are several 
cases (i.e., with very low optical luminosities and/or un- 
usual X-ray/optical flux ratios) that may be isolated X- 



ray binaries. If the counterparts are in fact (bound) clus- 
ters, then the observed X-ray luminosities arc well within 
the Eddington luminosities implied by the predicted cen- 
tral black hole masses (i.e., ~ 0.1% of the cluster mass 
estimated from the optical luminosity). We also find that 
the scenario in which the optical fiux is dominated by 
the emission of an accretion disk around a black hole re- 
sults in optical and X-ray fluxes consistent with those 
observed, although the black hole mass would have to be 

> IO'^Mq unless the disk is irradiated. 
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Appendix 

8. NOTES ON INDIVIDUAL ULXS 

8.1. ULX 1 

ULX 1 is found in disk of NGC 278 , which is a face-on 
spiral (T=3). 6 potential counterparts are fotmd. mostly 
in the 1-2ct error annulus although the brightest counter- 
part is within the la error circle. The B-V and V-I colors 
are generally blue although some counterparts have B-V 

> 0.5. 

8.2. ULX 2 

ULX 2 is found in NGC 891, an edge-on spiral galaxy 
(T=3), located just outside of a prominent dust lane. 
B-V = -0.2 and V-I < O.l, suggesting a very blue coun- 
terpart or else the optical source varied between the dif- 
ferent filter observations. 

8.3. ULXs 3-4 

These ULXs arc found in NGC 1068, a well-known 
face-on spiral galaxy (T=:3) with strong starburst and 
Seyfcrt activity. In both cases the source density is high 
and in most cases different counterparts arc found in dif- 
ferent filters. However in this case the astrometry could 
not be corrected, due in part to the complexity of the nu- 
clear region in both the X-ray and optical bands and the 
small FOV of the WFPC2 images. Future work may im- 
prove the registration by leveraging ground-based optical 
images. 

8.4. ULXS 

ULX is found in the disk of NGC 1097, a face-on spiral 
galaxy (T=3). The only Association image was in the 
UV and no counterparts were found. 

8.5. ULXs 6-8 

These ULXs are found in NGC 1399, an elliptical 

galaxy (T=-5). ULXs 7 and 8 have single counterparts 
that are consistent with globular clusters (see also An- 
gelini et al. 2001). ULX 6 is possibly due to a background 
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AGN, although its log Lx/Lv,eff value of >1.4 is some- 
what large relative to typical Chandra Deep Field sources 
(see Bauer et al. 2004). 

8.6. ULXs 9-11 

These ULXs are found in NGC 1700, an elliptical 
galaxy (T=-5). NGC 1700 is one of the most distant 
galaxies in our sample at 52 Mpc. No counterparts are 
found for ULX 10. The counterparts for ULX 9 and 11 
have red colors, consistent with a globular cluster inter- 
pretation. 

8.7. ULXs 12-14 

These ULXs are found in M82 (T=90), which is the 
archetypal, and one of the nearest, starburst galaxies. 
In all cases the optical emission is too complex for indi- 
vidual sources to be identified, with the source regions 
dominated by diffuse emission with non-uniform extinc- 
tion. 

8.8. ULXs 15-16 

ULXs 15 and 16 arc fomid in NGC 3256, a merger 
(T=99). Two counterparts are found for each ULX. One 
of the counterparts for ULX 16 has a red color limit (B-I 
> 1.0). This may be due to extinction however note that 
this counterpart is only marginally within the 2a error 
circle and may be unrelated to the ULX. 

8.9. ULXs 17-18 

These ULXs are in NGC 3311, the central galaxy of 
the Hydra cluster (T=-4) . The ULXs are found near the 
nucleus which is dusty. Multiple counterparts are found 
for both, with one being very red, most likely due to dust 
extinction. 

8.10. ULX 19 

ULX 19 is found in the edge-on spiral galaxy NGC 
3628 (T=3). The ULX is in the center of a promi- 
nent dust lane, and a single counterpart is detected with 
log Lx/Lv^ff = 1.6. This is consistent with an isolated 
X-ray lainary scenario however the V-band flux is proba- 
bly underestimated due to extinction. 

8.11. ULX 20 

ULX 20 is found in the spiral galaxy (T=3) NGC 3627, 
located on the edge of the disk. Here it should be noted 
that we were not able to correct the HST astrometry and 
the position of the ULX (relative to nearby stars) differs 
by ~ 6" between the two V band images. Therefore 
the two 3.3" {2a) error circles do not overlap and the 
detection of a counterpart (beyond the la error radius) 
in one image but not the other does not imply variability. 

8.12. ULX 21-28 

These ULXs are found in the Antennae galaxy merger. 
As discussed in Section 1, the ULX/HST correlation in 

the Antennae has been assessed in Zezas et al. (2002) 
and Zezas & Fabbiano (2002). They did not correct the 
HST astrometry and assumed a 2" error circle, and of- 
ten found multiple potential counterparts, typically with 
blue colors suggesting a young population. Our results 
arc consistent, with our (Method 2) 2a error circle being 
1.6" in radius. ULX 25 was identified as a background 
AGN by Clark et al. (2005). 



8.13. ULX 29 

ULX 29 is found in NGC 4150, a lenticular galaxy (T=- 
2). The only Association image available was an I band 
exposure in which no counterparts were found. 

8.14. ULX 30 

ULX 30 is found in NGC 4490 which is a late- type spi- 
ral galaxy (T=7) with the ULX being located along the 
disk. Not surprisingly, multiple HST sources are found 
within the 2a error circle. 

8.15. ULX 31 

ULX 31 is found in NGC 4559, also a late-type spiral 
(T=6). The ULX is located in the disk near the nucleus. 
Two faint point-like counterparts are found although the 
region is somewhat complex. Cropper et al. (2004) found 
no counterparts in the vicinity of this ULX (their XIO), 
using Association data. Cropper et al. (2004) also use 
B and I band images which were not available at the 
start of this study. Similarly the HST data discussed in 
Cropper et al. (2004) and Soria et al. (2005) overlapping 
another ULX in NGC 4559 (their X7) was not available 
at the start of this study. 

8.16. ULX 32 

ULX 32 is found in NGC 4565, an edge-on spiral 
(T=3). The ULX is on the edge of the bulge and has been 
identified with a globular cluster in Wu et al. (2002). 

8.17. ULX 33 

ULX 33 is found in NGC 4594 (the Sombrero galaxy; 
T=l). The ULX is located in the disk near the nucleus 
(outside the dust lane). An Ha counterpart is found 
marginally within the 2a error circle (a cosmic ray not 
completely removed is noticeable in the la error circle). 

8.18. ULX 34 

ULX 34 is found in NGC 4725, an early-type spiral 
(T=2). The ULX is located inside an outer spiral arm. 
Diffuse flux is evident at the ULX position and a (point- 
like) V-band counterpart is found at the 2a radius. 

8.19. ULX 35 

ULX 35 is fomid in NGC 4945, a late-type spiral 
(T=6). No counterpart is found, however the ULX re- 
gion is located near the nucleus within the disk and the 
lack of a counterpart may be due to extinction. 

8.20. ULXs 36-37 

These ULXs are found in the elliptical galaxy NGC 
5018 (T=-5), with ULX 37 being located near the nu- 
cleiis. The counterpart for ULX 37 has a very blue color 
(U-V=-1.6). This counterpart is therefore likely to be 
a background source, with an X-ray/optical flux ratio 
consistent with AGN. 

8.21. ULX 38-39 

These ULXs are found in the face-on, late-type spiral 
galaxy NGC 5033 (T=5). The ULXs are located close 
to each other, within the disk. Only UV Association 
images were available and several potential counterparts 
were detected. 
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8.22. ULX 40 

ULX 40 is found in the elliptical galaxy NGC 5128 
(Centaurus A; T=-2). The ULX is located outside of 
the dusty central region. Two counterparts are found. 
Neither of our counterparts corresponds to the counter- 
parts discussed in Ghosh et al. (2006), based on the same 
WFPC2 data. This is because Ghosh et al. applied as- 
trometric corrections to the HST data based on a cata- 
loged globular cluster (an approach not attempted here) 
and a nearby transient source that was not detected in 
the ACIS image used here (but was detected in a sub- 
sequent ACIS observation not public at the start of our 
study). The brighter of the two counterparts Usted here 
is ~ 1 magnitude brighter than the counterpart discussed 
in Ghosh et al., however we derive a similar color (V-I 
~ 1 compared with V-I ^ 0.7 in Ghosh et al.) implying 
this counterpart is similar in nature to the counterpart 
discussed in Ghosh et al. Since our discussion is based on 
the properties of counterparts on average (and often the 
counterparts to a given ULX have similar colors), our 
conclusions are not affected our (potentially incorrect) 
counterpart selection here. 

8.23. ULXs 41-42 

These ULXs are found in the Circinus galaxy, a face- 
on spiral (T=3). The ULXs are located in the disk of 
the galaxy. Only narrow-band Association images were 
available that overlapped ULX 41. ULX 41 is CG X-2 
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image was in the UV/U (F300W) and a single counter- 
part was detected. 
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TABLE 2 
HST/WFPC2 Datasets 

Note. — The tabic lists basic properties of the HST/WFPC2 data used for our study, selected as described in Section 2.2. Column (1) lists 
the HST/WFPC2 association name. Column (2) lists the host galaxy name. Column (3) lists the filter with which the data were obtained. When 
multiple datasets with the same filter are available for a given galaxj-. a unique identifier is added in parentheses. Column (4) lists the number of 
exposures in the WFPC2 association. Column (5) lists the total exposure time. Column (6) indicates the relative exposure times of the different 
exposures in the association. The first number is the shortest exposure time divided by the average exposure time. The second number is the 
longest exosure time divided by the average. Column (7) indicates the method by which the astrometry was calibrated, as described in Section 3. 
Method 1 coordinates were shifted to align the galaxy nucleus with the AGN coordinates in Chandra data; method 2 coordinates were shifted to 
align observed USNO-Bl.O stars with their catalog coordinates; and method 3 coordinates are obtained from the original HST pipeline headers. 
When the number is followed by the name of a dataset in square brackets, that means that the image was aligned with the listed dataset for the 
same galaxy (which itself was aligned using the listed method). For method 2 the number of USNO-Bl.O stars that was identified is listed in 
parentheses. The relative accuracy of the Chandra/HST alignment for each method is discussed in Section 3. Column (8) lists the ID numbers 
(from Table 1) of the ULXs that fall on the image. Columns (9)-(ll) list the (X,Y) pixel position in WFPC2 detector coordinates and WFPC2 
chip number that correspond to the position of the ULX given in Table 1. 
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F658N 


M 


23.34 


0.14 


-0.11 


0.18 


19.54 


0.03 


-12.14 


129.0 


0.3/29.4 


2 


F336W 


M 


24.69 


-0.49 


0.31 


0.58 


19.86 


0.02 


-11.82 


335.9 


0.3/32.6 




F555W 


M 


23.93 


-0.28 


0.23 


0.36 


20.37 


0.05 


-11.31 


96.0 


0.2/24.8 




F658N 


M 


23.34 


-0.17 


0.53 


0.55 


19.69 


0.04 


-11.99 


111.4 


0.2/24.8 


3 


F336W 


M 


24.69 


1.11 


0.94 


1.46 


21.22 


0.04 


-10.46 


95.1 


0.3/35.1 




F555W 


M 


23.92 


1.28 


0.84 


1.53 


21.30 


0.08 


-10.38 


39.3 


0.3/29.5 




F658N 


M 


23.34 


0.85 


0.94 


1.26 


19.64 


0.04 


-12.04 


116.6 


0.5/55.1 


4 


F336W 


M 


24.69 


-0.59 


0.15 


0.61 


20.78 


0.03 


-10.90 


142.7 


0.3/32.6 




F555W 


N 


21.81 
















1 


F336W 


M 


24.03 


-0.16 


-0.31 


0.35 


23.39 


0.25 


-8.29 


6.8 


0.2/15.8 




F555W 


N 


23.66 
















2 


F336W 


M 


24.03 


-0.26 


0.22 


0.34 


23.60 


0.30 


-8.08 


5.5 


0.2/18.9 




F555W 


N 


23.66 
















3 


F336W 


M 


24.03 


0.89 


-0.14 


0.90 


23.72 


0.34 


-7.96 


5.0 


0.2/21.0 




F555W 


N 


23.66 
















4 


F336W 


M 


24.03 


0.96 


0.78 


1.24 


23.44 


0.26 


-8.24 


6.5 


0.2/21.0 




F555W 


N 


23.66 
















5 


F336W 


M 


24.03 


0.67 


-1.32 


1.48 


23.81 


0.36 


-7.87 


4.5 


0.2/17.9 




F555W 


N 


23.66 
















6 


F555W 


N 


23.66 




















F658N 


M 


23.35 


-0.52 


0.40 


0.65 


21.84 


0.15 


-9.84 


15.0 


0.2/22.3 


1 


F336W 


M 


24.43 


1.19 


-0.77 


1.42 


21.49 


0.05 


-10.19 


57.4 


0.2/18.9 




F555W 


N 


23.51 


















F658N 


N 


22.77 
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TABLE 3 

Point-like Optical Sources Coincident with ULXs 



24 


NGC 4038 


1 


F336W 


M 


24.35 


0.84 


1.12 


1.40 


18.88 


0.01 


-12.80 


611.6 


0.2/26.2 








F555W 


c 






















2 


F658N 


M 


22.44 


1.04 


0.25 


1.07 


20.15 


0.07 


-11.53 


32.7 


0.3/27.5 


25 


NGC 4038 


1 


F336W 


M 


24.67 


-0.39 


0.17 


0.43 


20.58 


0.04 


-11.10 


157.0 


0.3/27.4 








F658N 


M 


23.28 


-0.07 


0.46 


0.47 


20.42 


0.07 


-11.26 


49.2 


0.2/24.4 


26 


NGC 4038 




F336W 
F555W 
F658N 


c 
c 
c 


















27 


NGC 4038 




F336W 
F555W 
F658N 


N 
N 
N 


24.62 
23.47 
23.44 


















28 


NGC 4038 




F555W 


N 


23.47 


















29 


NGC 4150 




F814W 


N 


23.85 


















30 


NGC 4490 


1 


F606W 


M 


24.44 


-0.57 


-1.00 


1.15 


24.15 


0.34 


-5.31 


4.9 


0.2/9.5 






2 


F606W 


M 


24.44 


0.98 


-0.12 


0.99 


23.08 


0.14 


-6.38 


13.4 


0.2/8.2 






3 


F606W 


M 


24.44 


0.56 


-1.11 


1.24 


23.52 


0.20 


-5.94 


8.9 


0.2/7.0 






4 


F606W 


M 


24.44 


0.02 


-1.50 


1.50 


23.76 


0.24 


-5.70 


7.1 


0.2/5.7 






5 


F606W 


M 


24.44 


-0.23 


1.28 


1.30 


23.54 


0.20 


-5.92 


8.7 


0.2/6.3 






6 


F606W 


M 


24.44 


1.11 


-1.30 


1.71 


21.40 


0.04 


-8.06 


64.3 


0.2/7.3 


31 


NGC 4559 


1 


F606W 


M 


24.50 


-0.62 


0.19 


0.65 


24.39 


0.37 


-5.54 


4.4 


0.2/8.6 






2 


F606W 


M 


24.51 


0.88 


-0.64 


1.09 


23.69 


0.21 


-6.24 


8.4 


0.2/10.2 


32 


NGC 4565 


1 


F450W(1) 


s 


25.64 


0.20 


0.12 


0.23 


23.81 


0.12 


-7.26 


20.7 


0.2/18.3 








F450W(2) 


N 


25.22 






















F555W 


N 


25.22 
























F814W(1) 


s 


24.78 


0.19 


0.13 


0.23 


22.55 


0.09 


-8.52 


30.7 


0.2/18.2 








F814W(2) 


N 


24.25 
















33 


NGC 4594 


1 


F658N 


s 


22.07 


0.53 


0.29 


0.60 


19.09 


0.04 


-11.72 


58.7 


0.3/18.0 


34 


NGC 4725 


1 


F606W 


M 


25.41 


0.15 


-0.59 


0.61 


22.08 


0.03 


-8.95 


76.2 


0.2/16.9 


35 


NGC 4945 




F606W 


N 


24.17 
















36 


NGC 5018 




F336W 
F555W 


N 
N 


23.30 
26.42 


















37 


NGC 5018 


1 


F336W 


M 


23.77 


0.01 


-0.05 


0.05 


22.00 


0.11 


-10.86 


17.9 


0.2/35.1 








F555W 


M 


25.67 


-0.05 


-0.02 


0.05 


23.58 


0.09 


-9.28 


22.8 


0.2/33.1 


38 


NGC 5033 


1 


F218W 


N 


21.86 






















F300W 


S 


24.18 


-0.22 


-0.07 


0.24 


23.91 


0.34 


-7.45 


5.1 


0.2/18.3 


39 


NGC 5033 


1 


F300W 


s 


24.19 


-0.23 


0.07 


0.24 


21.91 


0.08 


-9.45 


32.5 


0.2/19.5 






2 


F218W 


S 


22.06 


0.17 


-0.19 


0.26 


20.65 


0.14 


-10.71 


14.4 


0.2/17.2 


40 


NGC 5128 


1 


F555W 


s 


24.40 


-0.15 


-0.87 


0.88 


22.89 


0.18 


-5.56 


11.0 


0.3/6.4 








F814W 


S 


23.26 


-0.16 


-0.85 


0.87 


21.99 


0.21 


-6.46 


8.7 


0.2/5.7 






2 


F555W 


N 


24.32 






















F814W 


M 


23.26 


0.73 


-1.04 


1.27 


22.28 


0.27 


-6.17 


6.5 


0.2/5.0 


41 


Circinus 


1 


F502N 


M 


22.67 


0.50 


0.17 


0.53 


20.70 


0.10 


-7.14 


23.5 


0.2/3.5 








F656N 


M 


21.92 


0.42 


0.25 


0.49 


17.93 


0.03 


-9.91 


154.5 


0.2/3.2 


42 


Circinus 




F502N 
F606W 
F656N 


N 
N 
N 


22.30 
24.97 
21.67 
















43 


NGC 6500 




F218W 
F547M 


N 
X 


21.53 


















44 


NGC 7174 


1 


F300W 


M 


23.28 


0.39 


0.14 


0.42 


22.54 


0.25 


-10.30 


6.6 


0.2/42.7 



Note. — Column (1) lists the ULX ID number from Table 1. Column (2) lists the host galaxy name. Column (3) lists a running ID number of 
the optical sources spatially coincident with the ULX given in column (1). These optical source ID numbers are also used to label the sources in 
Figure 2. Column (4) lists the filter used for the HST observation, as in Table 2: the corresponding HST dataset name can be obtained from that 
table. Column (5) indicates the optical morphology of the region surrounding the ULX, as defined in Section 5. Column (6) lists the 4cr limit for 
source detection determined as described in Section 5. Columns (7) and (8) lists the offset in RA and Dec. of the optical sources to the ULX in 
arcseconds. Column (9) lists the total offset of the optical sources to the ULX in arcseconds. Column (10) lists the apparent magnitude of the 
optical source in the VEGAMAG system. Column (11) lists the formal error in the apparent magnitude. Column (12) lists the absolute magnitude 
of source, calculated using the galaxy distance in Table 1. Column (13) lists the S/N level of the optical source detect ion. Column (14) lists an 
estimate of the source extent in arcseconds and pc assuming the distances listed in Table 1 (in most cases the sources are unresolved in which case 
this estimate is intended to give an upper-limit to the physical source size). 



TABLE 4 
WFPC2 Counterpart Colors 



ULX 


Galaxy 


T-type 


Src. 


Color 


Value 


ID 




RC3 


No. 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


1 


n0278 


3.0 


1 


B - V606 


0.75 ±0.12 


1 


n0278 


3.0 


1 


B - I 


1.40 ±0.12 


1 


n0278 


3.0 


1 


Veoe - I 


0.65 ± 0.09 


1 


n0278 


3.0 


2 


B - Veoe 


0.50 ±0.22 


1 


n0278 


3.0 


2 


B - I 


1.26 ±0.20 


1 


n0278 


3.0 


2 


Veoe - I 


0.76 ±0.18 


1 


n0278 


3.0 


3 


B - Veoe 


0.28 ± 0.08 


1 


n0278 


3.0 


3 


B - I 


0.66 ± 0.09 


1 


n0278 


3.0 


3 


Veoe - I 


0.38 ± 0.09 
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TABLE 4 — 


Continued 




ULX 


1-1 a 1 n v\7 

\jcxi.<XiA.y 




Src. 


Color 


V 1 


ID 




RC3 


No. 


















1 


n0278 


3.0 


4 


B - V(jo6 


-0.01 ± 0.44 


1 


n0278 


3.0 


4 


B - I 


0.33 ±0.48 


1 


n0278 


3.0 


4 


Vgog - I 


0.34 ±0.54 


1 


n0278 


3.0 


5 


B - V(i()fi 


0.42 ±0.15 


1 


n0278 


3.0 


5 


B - I 


0.95 ± 0.15 


1 


n0278 


3.0 


5 


Veoe " I 


0.53 ±0.14 


1 


n0278 


3.0 


6 


Veoe - I 


0.65 ± 0.24 


1 


n0278 


3.0 


6 


B - Veoe 


> 1.52 


1 


n0278 


3.0 


6 


B - I 


> 2.17 


2 


n0891 


3.0 


1 


Veoe - I 


< 0.13 


2 


n0891 


3.0 


1 


B - Veoe 


-0.17 ±0.19 


2 


n0891 


3.0 


1 


B - I 


< -0.04 


7 


nl399 


-5.0 


1 


B - I 


1.65 ±0.04 


7 


nl399 


-5.0 


1 


B - Veoe 


0.53 ± 0.04 


7 


nl399 


-5.0 


1 


Veoe - 1 


1.12 ±0.03 


8 


nl399 


-5.0 


1 


B - I 


1.87 ±0.03 


8 


nl399 


-5.0 


1 


B - Veoe 


0.91 ±0.03 


8 


nl399 


-5.0 


1 


Veoe - I 


0.96 ± 0.02 


9 


nl700 


-5.0 


1 


V555 - I 


3.03 ± 0.03 


11 


nl700 


-5.0 


1 


V555 - I 


1.80 ±0.08 


11 


nl700 


-5.0 


1 


V555 - Veoe 


0.50 ±0.07 


11 


nl700 


-5.0 


1 


Veoe - I 


1.30 ±0.04 


15 


n3256 


99.0 


1 


B - I 


0.79 ±0.28 


15 


n3256 


99.0 


1 


B - V555 


0.37 ± 0.36 


15 


n3256 


99.0 


1 


V666 - I 


0.42 ± 0.37 


16 


n3256 


99.0 


2 


B - I 


> 1.03 


17 


n3311 


-4.0 


1 


V666 - I 


> 2.45 


17 


n3311 


-4.0 


3 


V666 - I 


1.12 ±0.18 


17 


n3311 


-4.0 


4 


V556 - I 


< 0.08 


18 


n3311 


-4.0 


1 


V555 - I 


< 0.66 


18 


n3311 


-4.0 


2 


V555 - I 


1.10 ±0.11 


21 


n4038/9 


9.0 


1 


U - V555 


< 0.34 


21 


n4038/9 


9.0 


2 


U - V555 


-0.51 ± 0.05 


21 


n4038/9 


9.0 


3 


U - V565 


-0.08 ± 0.09 


21 


n4038/9 


9.0 


4 


U - V556 


< -1.03 


22 


n4038/9 


9.0 


1 


U - V555 


< -0.27 


22 


n4038/9 


9.0 


2 


U - V555 


< -0.06 


22 


n4038/9 


9.0 


3 


U - V555 


< 0.06 


22 


n4038/9 


9.0 


4 


U - V555 


< -0.22 


22 


n4038/9 


9.0 


5 


U - V555 


< 0.15 


23 


114038/9 


9.0 


1 


U - V555 


< -2.02 


32 


n4565 


3.0 


1 


V555 - 1 


> 2.67 


32 


n4565 


3.0 


1 


B - I 


1.26 ±0.15 


32 


n4565 


3.0 


1 


B - VsBB 


< -1.41 


37 


n5018 


-5.0 


1 


U - V665 


-1.58 ±0.14 


40 


n5128 


-2.0 


1 


V555 - 1 


0.90 ±0.27 


10 


115128 


-2.0 


2 




2.01 



Note. — Column (1) lists the ULX ID number from Tabic 1. Column (2) lists the host galaxy name. Coh.imn (3) lists the galaxy morphology 
T-type from the Third Reference Catalog of Bright Galaxies (RC3, de Vaucouleurs et al. 1991). Column (4) lists the id number of the source from 
Table 3. Column (5) lists the color being computed, with the filters designations U, B, V555, Veoe and I refer to HST WFPC2 filters F336W, 
F450W, F555W, F606W and F814W. Column (6) lists the color or color limit. 



TABLE 5 

WFPC2 Counterpart Optical Luminosities 



ULX 


Src. 


log Lx/Lq^X 


Veff 


Mv,eff 


log Lv,eff/LQ,V 


logLx/Lv,eff 


ID 


No. 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


1 


1 


5.4(B) , 5.4(V) , 5.4(1) 


21.9 


-8.5 


5.3 


0.5 


1 


2 


5.1(B) , 5.1(V) , 5.1(1) 


22.7 


-7.6 


5.0 


0.9 


1 


3 


5.7(B) , 5.5(V) , 5.4(1) 


21.3 


-9.0 


5.5 


0.3 


1 


4 


4.9(B) , 4.7(V) , 4.5(1) 


23.4 


-6.9 


4.7 


1.2 


1 


5 


5.3(B) , 5.3(V) , 5.2(1) 


22.0 


-8.4 


5.3 


0.6 


1 


6 


5.0(V) , 4.9(1) , <4.6 (B) 
4.6(V) , 4.9(B) , <4.4 (I) 


23.0 


-7.4 


4.9 


1.0 


2 


1 


23.2 


-6.8 


4.6 


1.2 


6 


1 


<4.5 (I) , <4.4 (B) , <4.4 (V) 
5.7(1) , 5.6(B) , 5.6(V) 


>25.9 


>-5.5 


<4.1 


>1.4 


7 


1 


22.6 


-8.8 


5.5 


0.2 


8 


1 


6.0(1) , 5.8(B) , 5.9(V) 


21.8 


-9.6 


5.8 


-0.0 


9 


1 


6.5(V) , 7.4(1) 


22.1 


-11.5 


6.5 


-0.6 


10 


1 


<5.2 (V) , <5.4 (I) , <5.4 (V) , <5.8 (I) 


>25.4 


>-8.1 


<5.2 


>0.8 


11 


1 


5.8(V) , 6.2(1) , 6.0(V) , 6.0(V) 


23.6 


-10.0 


5.9 


0.4 


15 


1 


5.3(B) , 5.1(1) , 5.1(B) , 5.2(V) 


24.6 


-8.2 


5.2 


0.4 
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ULX Src. 
ID No. 

(1) (2) 



(3) 



Veff Mv,eff ^OgLv,eff/LQ,v ^OgLx/Lv,eff 



(4) 



(5) 



(6) 



(7) 



15 


2 


5.2(1) , <4.9 (I) , <5.1 (I) , <5.1 (V) , <5.1 (I) 


24.1 


-8.8 


5.4 


0.2 


16 


1 


<5.1 (B) , <4.9 (I) 


>25.6 


>-7.2 


<4.8 


>1.1 


16 


2 


5.0(1) , <5.1 (B) 


25.1 


-7.7 


5.0 


0.9 


17 


1 


5.5(1) , <4.8 (V) 


24.5 


-9.0 


5.5 


0.7 


17 


2 


5.4(V) 


24.8 


-8.7 


5.4 


0.8 


17 


3 


5.3(V) , 5.4(1) 


25.1 


-8.4 


5.3 


0.9 


17 


4 


5.3(V) , <5.0 (I) 


25.1 


-8.4 


5.3 


0.9 


18 


1 


5.1(V) , <5.0 (I) 


25.7 


-7.8 


5.1 


1.5 


18 


2 


5.6(V) , 5.8(1) 


24.2 


-9.3 


5.6 


0.9 


19 


1 


4.2(V) 


23.7 


-5.8 


4.2 


1.6 


20 


1 


4.1(V) , <3.6 (V) , <3.4 (I) 


23.6 


-5.5 


4.1 


1.5 


21 


1 


6.0(U) , <5.9 (V) 


22.2 


-9.4 


5.7 


0.1 


21 


2 


6.9(U) , 6.4(V) 


20.4 


-11.3 


6.4 


-0.6 


21 


3 


6.4(U) , 6.1(V) 


21.3 


-10.4 


6.1 


-0.2 


21 


4 


6.6(U) , <5.9 (V) 


20.9 


-10.8 


6.2 


-0.4 


22 


1 


5.5(U) , <5.1 (V) 


23.5 


-8.2 


5.2 


1.1 


22 


2 


5.4(U) , <5.1 (V) 


23.7 


-8.0 


5.1 


1.2 


22 


3 


5.4(U) , <5.1 (V) 


23.8 


-7.9 


5.1 


1.2 


22 


4 


5.5(U) , <5.1 (V) 


23.5 


-8.2 


5.2 


1.1 


22 


5 


5.4(U) , <5.1 (V) 


23.9 


-7.8 


5.0 


1.2 


22 


6 


<5.1 (V) 


>23.7 


>-8.0 


<5.1 


>1.2 


23 


1 


6.3(U) , <5.2 (V) 
7.3(U) 


21.6 


-10.1 


6.0 


-0.4 


24 


1 


19.0 


-12.7 


7.0 


-1.5 


25 


1 


6.6(U) 


20.7 


-11.0 


6.3 


-0.3 


27 


1 


<5.0 (U) , <5.2 (V) 


>24.7 


>-7.0 


<4.7 


>1.4 


28 


1 


<5.2 (V) 


>23.5 


>-8.2 


<5.2 


>0.3 


29 


1 


<4.0 (I) 


>24.0 


>-6.0 


<4.3 


>1.3 


30 


1 


4.0(V) 


24.1 


-5.3 


4.0 


1.5 


30 


2 


4.5(V) 


23.1 


-6.4 


4.5 


1.1 


30 


3 


4.3(V) 


23.5 


-5.9 


4.3 


1.3 


30 


4 


4.2(V) 


23.8 


-5.7 


4.2 


1.4 


30 


5 


4.3(V) 


23.5 


-5.9 


4.3 


1.3 


30 


6 


5.1(V) 


21.4 


-8.1 


5.1 


0.4 


31 


1 


4.1(V) 


24.4 


-5.5 


4.1 


1.9 


31 


2 


4.4(V) 


23.7 


-6.2 


4.4 


1.6 


32 


1 


5.0(1) , 5.1(B) , <4.3 (V) , <4.3 (I) , <4.5 (B) 
5.5(V) 


23.3 


-7.7 


5.0 


1.5 


34 


1 


22.4 


-8.7 


5.4 


0.4 


35 


1 


<3.7 (V) 


>24.2 


>-4.4 


<3.7 


>1.9 


36 


1 


<6.0 (U) , <4.5 (V) 


>26.5 


>-6.4 


<4.5 


>1.8 


37 


1 


6.5(U) , 5.6(V) 


23.6 


-9.3 


5.6 


0.3 


40 


1 


4.2(1) , 4.1(V) 


22.9 


-5.6 


4.1 


1.4 


40 


2 


4.1(1) , <3.6 (V) 


22.9 


-5.5 


4.1 


1.4 


42 


1 


<3.1 (V) 


>25.0 


>-2.9 


<3.1 


>3.0 



Note. — Column (1) lists the ULX ID number from Table 1. Column (2) lists ID number of the optical source from Table 3. Column (3) lists 
the optical luminosity computed for the band(s) in which photometric data was available. Column (4) lists Vg//, the effective V band magnitude 
computed based on the observed or assumed color of the source. Column (5) lists the absolute magnitude computed using Veff and the distance 
listed in Table t:sample. Column (6) lists log Lv,e///-^0,V) the optical luminosity of the source based on the Veff estimate, in solar luminosities. 
Column (7) lists the ratio of X-ray to optical luminosity. 



